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A quantum particle can be localized in a disordered potential, the effect known as Anderson local- 
ization. In such a system, correlations of wave functions at very close energies may be described, due 
to Mott, in terms of a hybridization of localized states. We revisit this hybridization description and 
show that it may be used to obtain quantitatively exact expressions for some asymptotic features of 
correlation functions, if the tails of the wave functions and the hybridization matrix elements are as- 
sumed to have log-normal distributions typical for localization effects. Specifically, we consider three 
types of one-dimensional systems: a strictly one-dimensional wire and two quasi-one-dimensional 
wires with unitary and orthogonal symmetries. In each of these models, we consider two types of 
correlation functions: the correlations of the density of states at close energies and the dynamic 
response function at low frequencies. For each of those correlation functions, within our method, we 
calculate three asymptotic features: the behavior at the logarithmically large "Mott length scale", 
the low-frequency limit at length scale between the localization length and the Mott length scale, 
and the leading correction in frequency to this limit. In the several cases, where exact results are 
available, our method reproduces them within the precision of the orders in frequency considered. 

PACS numbers: 73.20.Fz, 73.21.Hb, 73.22.Dj 



I. INTRODUCTION 

The localization of a quantum particle in a disordered 
potential (commonly known as Anderson localizatiorP^ is 
one of the most fascinating mesoscopic phenomena (see, 
e.g., Refs. |2] and [3| for a review). Arising from quan- 
tum interference between different particle trajectories, 
localization depends strongly on the dimensionality of 
the system. In one dimension, such an interference is 
most relevant, and an arbitrarily weak potential is known 
to localize a particle (in the absence of decoherence)Pt^. 
Besides, one-dimensional case is most accessible for ana- 
lytic studies, which makes it the best understood model 
of localization (see, e.g., Refs. [5] and IH]) . 

For the purpose of the present paper, we distin- 
guish several models of one-dimensional localization: the 
strictly-one-dimensional (SID) case (with one conduct- 
ing channel) and the quasi- one- dimensional (QID) wire 
(with ^ I conducting channels). These two limits 
exhibit some common universal properties, but are typi- 
cally treated with different analytic techniques (Berezin- 
sky techniqu d^" * ^^* and contemporary method s-^ -'^ in the 
SID case, and the sigma- model techniqu^^^ti^ in the 
QID case). The quasi-one-dimensional wires may be fur- 
ther classified in terms of the symmetries of the Hamil- 
tonian, according to the rand om-matrix-theory scheme 
(unitary, orthogonal, etc.)P^E3. 

One of the main quantitative characteristics of the lo- 
calization is the statistics of (localized) eigenfunctions. 
In one dimension, it was stu died extensively, and many 
analytic results are 

availabltPEsHUl. Most of the ana- 
lytical results are derived in the weak-disorder regime 
(which is believed to obey the single-parameter-scaling 
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FIG. 1: A schematic (not numerically exact) view of the cor- 
relation functions R{u),r) (top panel) and S{uj,r) (bottom 
panel) defined in Eqs. ([T]) and H, respectively. The dashed 
lines denote the oj — > limits [the same function (361 for 
R[lo, r) and S{ijj, r)]. The shaded regions denote SR{lu, r) and 
5S{ijj, r) as defined in Eq. (j4|. The features at the Mott length 
scale Lm are denoted by Ruii^, r) and Sm{lj, r), respectively. 



propert; 



see also Refs. '25', 1261 and [8] for further dis- 



cussions). Remarkably, the statistics of the "envelopes" 
of localized eigenfunctions in this regime is universal for 
SID and QID problems (independently of the symme- 
try class) and can be e xpre ssed in terms of the Liou- 
ville quantum mechanics^ES, while the short-range os- 
cillations distinguish between SID and QID cases and 
between symmetry classes in the QID case. 

A more detailed information about localization (in par- 
ticular, relevant to dynamic properties) can be extracted 
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TABLE I: Correlation function R(u},r) at low frequencies {uj <^ 1). 
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TABLE IL Correlation function S{u),r) at low frequencies (cj ^ 1). 
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SID 
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Q ID-orthogonal 
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from correlations between eigenfunctions at different en- 
ergies. Two such quantities may be defined^'"^^': the 
density-of-states (DOS) correlation function, 

Riuj, \xi - X2I) - ly-^ 6iEr, - E)6iE„, - E - lo) 

n,m 

X |^„(a;i)|>™(x2)p), (1) 
and the dynamic response function, 



Siuj, |xi - X2I) = i^-^ Q2 '^(^" - E)5{E^ ~E-u) 

X ■>Jj^{xi)Tjjn{x2W,ni^2)i^mixi)) . (2) 



Here the sum is taken over the eigenstates ipn with en- 
ergies En, and z/ is the average density of states. The 
normalization of these correlation functions is chosen in 
such a way that they are dimensionless quantities with a 
finite limit in an infinitely long wire. It will be further- 
more convenient to measure the lengths in the units of 
the localization length ^ and the energies in the units of 
the average level spacing within the localization length 
Ap^. With this convention, R{uj,r) and S{uj,r) become 
dimensionless functions of dimensionless parameters. 

Both r) and Sjuj , r) w ere studied analytically in 
detail in the SID modep32Sl^ and i?(aj,r) has been re- 
cently calculated in the QlD-unitary modeP^ (in all 
those studies, the weak-disorder limit was assumed). 
While the limiting form of these correlations at — > 
is determined by the single-wave-function statistics and 
is therefore universal for both SID and QID models, 
the corrections at finite uj distinguish between SID and 
QID''" 31^ Qualitatively, the properties of the correlation 
bmctions R(uj,r) and S{uj,r) in the low-frequency limit 
a; <C 1 may be understood using the original argument 
by Mott about the hybridization of the localized wave 
function^ (the opposite hmit lo ^ 1 can be studied by 
means of the standard perturbation theory). However, 
the first attempt to promote the Mott's arguments to 
quantitative calculations in Ref. produced some in- 
accurate results (as can be seen by comparing to exact 



expressionP^ , since it neglected mesoscopic fluctuations 
of the tails of the localized states. 

In the present work, we rectify this approach and re- 
visit Mott's arguments on wave-function hybridizatiorP^l 
taking into account the log-normal distribution of the 
tails of the localized states^. Our method is based on 
a number of assumptions that we introduce in the main 
text and then explicitly summarize and discuss in the last 



section of the paper (Section VIII I. As a result, we obtain 



a semi-phenomenological description of the hybridization 
of the localized states at distances much larger than the 
localization length, r 3> 1. Our theory reproduces cor- 
rectly the physics at the "Mott length scale" , 



Lm = 21n(l/a 



(3) 



and the leading correction to R{uj,r) at 1 ^ r < 
Lm in the QlD-unitary model. We further make pre- 
dictions concerning the properties of i?(a;,r) in the 
QlD-orthogonal model and of ^(a;, r) in all the above- 
mentioned models. These predictions may be checked 
against future sigma-model calculations in QID systems. 



II. MAIN RESULTS 

In the present work, we assume the single-parameter- 
scaling regime for the tails of the localized states. 
Namely, we suppose that at distances r ^ 1, the de- 
cay of the localized wave function may be described by a 
log-normal distribution with the width and median (or, 
formally, the variance and the mean of the logarithm) 
described by one parameter and with an appropriate cut- 
off of the tails. By combining this assumption with the 
Mott's argument about the wave-function hybridization 
(see subsequent sections for details), we can infer quanti- 
tative details about the behavior of the correlation func- 
tions ([ij and ([2| in the low-frequency limit a; ^ 1. 

The general structure of those two correlation func- 
tions contains two main separate regimes: r <C Lm and 
r ^ Lm (Fig. ^. At r <C Lm, the correlations are 
known to be dominated by the statistics of a single wave 
function-^^ and it is natural to represent them 
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as 



i?(w,r) 
5(w,r) 



R{u; 



0, r) 
0,r) 



-SR{Lu,r) 
6S{uj,r) , 



(4a) 
(4b) 



where 6R{Lj,r) and (5S'(aj,r) vanish as a; — > 0. 

At r ~ the correlation function R{uj,r) exhibits 
a crossover from zero to one centered at Lm and with 
a width of the order \/Lm, and the correlation fu nction 
S{uj,r) has a negative bump at the same locatiorp2i2fil. 
The asymptotic form of these features at w — will be 
denoted as Rm{^,t) and 5'm(w,''), respectively. 

Our hybridization argument reproduces the (univer- 
sal) main asymptotics R{uj — 0,r) and S{uj -> 0,r) at 
1 ^ r ^ Lm, the (nonuniversal) leading in w correc- 
tions SR{u!,r) and (5S(a;,r), as well as the universal be- 
havior of RM{uj,r) and Sm{'^,i'), see 

Tables |l] and |nf 4 

Some of these results can be verified against the exist- 
ing exact calculations, while others present new conjec- 
tures. As a byproduct of our calculation, we also relate 
the R{u) — >■ 0,r) and S{ijJ — 0,r) to the statistics of a 



single wave function [Eq. ( 39 ) below] , including the pro- 
portionality coefficient. 



III. STATISTICS OF WAVE-FUNCTION TAILS 

We start with a simplified statistical description of a 
single localized state in terms of the log-normal distri- 
bution of its tails. The statistics of a single wave func- 
tion has been studied in detail in both QID and SID 
geometried^EU^ and we first briefly summarize the exist- 
ing results and then propose our approximation. 

First of all, a localized state ^'(s^) can be represented 
as a product of a slowly varying envelope and a rapidly 
oscillating short-range coniponentpl 



(5) 



(here we include the dimensional factor (A^)^^/'^, where 
A is the cross section of the wire, in order to simplify 
the formulas below). The short-range component f{x) is 
correlated on the scale of the mean free path / and os- 
cillates on the scale of the particle wave length A^-. We 
choose it normalized to (|(/?(a;)p) = 1. The "envelope" 
component ip{x) is correlated on the scale of the local- 
ization length ^ (in SID, ^ ^ ^; in QID, £, ^ I) and does 

not oscillate. The two components ^'{x) and ^p{x) are dis- 
tributed independently. It was shown in Refs. I9I21I that 
such a decomposition is exact with the statistics of 
being universal for QID and SID systems, and that of 
(f{x) distinguishing between SID and QID and between 
different symmetry classes in QID. 

The statistics of ip{x) can be most conveniently de- 
scribed in terms of its logarithm 



As shown in Ref. [S] (section 3.2.2), the statistics oix{x) is 
given by a functional integral, which involves a diffusion- 
type quadratic action in xi^) a delta-function con- 
straint imposing the normalization of the wave function 
/ e^dx = 1. 

An accurate treatment of that functional integral is 
difficult, and we simplify it by observing that the tails 
of i/j^x) contribute very little to the normalization, and 
therefore the normalization-enforcing delta-function term 
is of little importance for the tails of ip{x). The normal- 
ization is mostly determined by the maxima of xi^), ^iid 
therefore the main role of this delta-term is to normalize 
the maxima of the function xi^)- We expect that the 
distribution of the maxima of xi^) ^ width of order 
one and centered around zero. 

This suggests our approximation for studying the 
wave-function tails. Instead of working with a full path 
integral of Ref.[9l we first fix the position xq and the value 
x{xq) of the maximum of x(a;) (with |x(a:o)| < 1) and 
replace the normalization constraint by an approximate 
condition that x(0) < x(3;o) everywhere. This guarantees 
a normalization of the wave function to a "logarithmic" 
precision: namely, the normalization of wave functions 
constructed in such a way will be of order one. Further- 
more, we will only be interested in a "coarse-grained" 
behavior of x(a;): typical scales of interest of xi^) will 
be of order r, and therefore for many purposes we do not 
need to distinguish between the maximum value xi^o) 
and zero. 

We thus arrive at the following coarse-grained descrip- 
tion of the ensemble of the envelopes x(^)' ^ localized 
state is determined by the location xq of its (global) max- 
imum (where xi^o) ~ 0) and the functional measure for 
the tails (which follows directly from the formula (3.34) 
of Ref. 9) 



exp 



dx 
dx 



+ sign{x ~ Xq) 



dx Dx{x) (7) 



with the constraint x{x) < for all x. In particular, the 
left and right parts of the tails {x < xq and x > xq, re- 
spectively) are distributed independently. The action ([T]) 
describes a diffusion with a drift, and the resulting form 
of the probability distribution for x{^) is approximately 
normal, with its variance growing linearly with x and its 
average decreasing linearly with x, sls x moves away from 
the center xq. 

For our calculations, we will be interested in one- and 
multi-point probability distributions of xi^) for ^ fixed 
position of xq. Let us start with the one-point probability 
distribution P(x,r), where r = \x — xo\. The action ([t]) 
results in the differential equation 



dP 
dr 



(8) 



X(x)=ln|7^(x)|2 



(6) 



d^P dP 
dx^ dx 

which describes diffusion with a drift. This equation, 
together with the boundary condition P{x > 0) = 0, 
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P(Z,r)^ 




FIG. 2: A description of a single-wave-function statistics. The 
wave function is described by the position a;o of its maximum 
and the distribution of its tails. The one-point probability 
distribution P{x, r) is assumed to have the form (|9| at r :§> 1. 



results in the long-time (r ^ 1) asymptotic form of the 
solution 

P{x,r)^f{^)Po{x,r), X<0, (9) 



where 



Poix,r) 



1 



exp 



4r 



(10) 



is the normal distribution (Fig. [2|. The effect of the 
boundary condition is the "cut-off" factor f{x/f)- The 
exact form of the function f{z) is determined by a short- 
time evolution (at r 1) and is therefore beyond our 
approximation scheme. The only property that we as- 
sume about /(z) (which becomes useful in Section VI I is 
its asymptotic behavior 



at 



-0. 



(11) 



The normalization of the probability distribution P(x, r) 
implies /(—I) = 1- Furthermore, the typical values of 
x{x) are not affected by the cut-off factor, and we have 
a "single-parameter-scaling" relation 



-((Ax)2)=-(X>='-. 



(12) 



Note that a similar single-parameter sc aling i s also well- 
known for the conductance of long wire^^^'^^JfUl 

The above consideration may be directly extended to 
multi-point probability distributions. For example, con- 
sider the probability distribution to find x(^i) = Xi and 
x{^2) — X2 under the condition that the maximum of 
x(a;) is located at xq (so that xi^o) ~ 0). The form of 
this probability distribution Pxoixii^ 11X2^X2) depends 
on the relative positions of xi, X2, and xq (see Fig. [s]). If 
Xi and X2 he on opposite sides of Xq (Fig. |3^), then the 
joint probability distribution factorizes: 



Pxo{Xi,Xi]X2,X2) = P{xi-,ri)P{x2-,r2) , 



(13) 



where — \xi — xo\ and P{x,r) is given by Eq. 

In other words, the left and right tails are statistically 

independent. In the opposite case, if xi and X2 belong 




Xi Xo X2 X 



Xq Xi X2 X 



FIG. 3: The two-point probability distribution 

Pxoixiy ^I'lX^j for ^ single wave function depends on the 
ordering of the observation points xi, X2, and the "center" of 
the wave function xq. (a) For xi and X2 belonging to differ- 
ent tails, the distribution factorizes, see Eq. (131. (b) For xi 



and X2 belonging to the same tail, the distribution is given 
by Eq. (fni. 



to the same tail, the distributions of x(a;i) and x{x2) are 
correlated. For the configuration shown in Fig. (the 
point xi lies between xg and 0:2), one gets 

Pxo{Xl,Xi;X2,X2) ^'(Xl,''l)^o(X2-Xl,?'2-?'l) ■ (14) 

Note that the second factor does not involve the cut-off 
function /(z), since, at — xi ^ 1 and — X2 ^ 1, the prob- 
ability of the functional integral ([t]) to return to xi^) — 
is exponentially small. 

We may note in passing that the two-point distribu- 
tions (131 and (14) are consistent with the one-point dis- 
tribution ([9]). Namely, 



dxi Pxo{xi,xi;x2,X2) = P{X2, \x2 - xq\) , (15) 



irrespectively of the relative positions of the points xq, 
xi, and X2- 

Generalization of this construction to many-point dis- 
tributions is straightforward. The only requirements are 
that the distances between all the points involved exceed 



the localization length, \xi — x-^ 
tails are considered, — Xi ^ 1. 



» 1, and that only small 



IV. WAVE-FUNCTION HYBRIDIZATION 

It was realized in the early works on localization that, 
at lli ^ Aj, the correlation functions ([T| and ([2| may be 
unde rstood in terms of the hybridization of two localized 
statep2133|_ Following the original argument, we may cut 
the wire into smaller pieces and consider two states ipA 
and ijjB localized in different pieces (centered at positions 
XA and xb, respectively, and with energies Ea and Eb)- 
As we connect the pieces of the wire together, the states 
get hybridized, and such pairs of states give the main 
contribution to the correlation functions ([I]) and ([2]). 

To use this argument at a quantitative level, we need 
to introduce the "hybridization" matrix element J be- 
tween the states ipA and ipB- Then the hybridized wave 
functions are given by the linear combinations 



^4 



-^B; (16) 
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FIG. 4: A schematic view of the hybridization of the two 
localized wave functions and ^b, as described by Eq. ( 16 1. 



where 



1 (^ ^ 



e = Eb — Ea and 



A = v^e^ +4|J|2 



(17) 



(18) 



are the energy splittings before and after hybridization 
(Fig. |4]). Such a pair of hybridized states contributes to 
the correlation functions ([T]) and ([2]), when A = w. 

It turns out that this approximate description repro- 
duces quantitatively many features of the exact results, 
provided the distribution of the tails ^ is taken into 
account, and appropriate assumptions on J are made. 

Namely, by analogy with the hybridization of states lo- 
calized in potential well^^, we assume that the hybridiza- 
tion matrix element J is proportional to the product of 
the two envelopes tpA{x) and iPb{x)'. 



\J\ = $VA(a;)V'B(a;), 



(19) 



where $ is a coefRcient of order one, which takes into 
account the short-range oscillations of the wave functions 
tpAix) and ipBix)- The distribution of $ is assumed to 
be statistically independent of the distributions of the 
envelopes iPa{x) and ^/'_b(x). The average ($) is taken to 
be of order one, so that Eq. ( 19 ) gives the matrix element 
J in the units of Aj. 

The specific properties of the distribution of $ will be 
of relevance for some of our calculations below. In fact. 



it is this distribution that distinguishes between the SID 
and QID geometries and between the different symmetry 
classes in the QID case. Specifically, in Section \VTI\ we 
will need the behavior of the probability distribution of <& 
at $ — >■ 0. Based on an analogy with the random-matrix 
theory, in that part of the calculation, we will use the 
following ansatz: 

dP($) =,5($-$o)d$ with$o--l in SID, (20a) 
dP{^) cx $ d$ , $ in QID unitary , (20b) 
dP{^) cx d$ , $ ^ in QID orthogonal . (20c) 

Our ansatz for (iP($) in the QlD-unitary and QID- 
orthogonal cases can be understood in terms of the sum 
of the hybridization amplitudes over a large number of 
channels. In the case of the unitary symmetry class, this 
sum is complex, and therefore its absolute value is dis- 
tributed as $ (i$ at small while in the orthogonal sym- 
metry class it is real with the measure at small 
The wave functions ipA^x) and iPb{x) in Eq. (19) are 



taken at some common point x in the tails of the wave 
functions. One can verify that, due to the log-normal 
statistics of the tails described in Section |III[ the prob- 
ability distribution of the product 7/1^1(2;) 7/1^(2;) is inde- 
pendent of the specific position of the point x. In other 



words, our ansatz ( 19 1 gives a consistent definition of the 
probability distribution of |J|. Equivalently, one may 
also rewrite 



|J| =$e'^'^/^ 



(21) 



where the parameter XJ has a distribution of the type 
with or without a cut-off factor (depending on the 

type of the points where the values of ipA (x) and ^pB (x) 
are fixed). 

We are now ready to formulate the improved version 
of the Mott hybridization argument by combining the 
three ingredients: (i) the hybridization of the wave func- 
tions ( [l6| , (ii) the statistical properties of a single wave 
function ([t]), and (iii) the properties of the hybridiza- 



tion matrix element (|19|). To obtain the uj <C Aj limits 



of the correlation functions ([T]) and ([2]), we restrict the 
sums over m and n to the two hybridized states ( 16 ) and 
arrive at 



i?(w,|xi -X2I) = / dxAdxB / dfi,^[xA{x)]dfi^JxBix)] / dP($) / de\ib+ixi)\^\^^ix2)\^SiA-u) , (22) 



S{iLi,\xi - X2\) = dxAdxB dfia:A[XA{x)]dfixB[XB{x)] dP{^) deip*^{xi)il;+{x2)ip*_(x2)ilJ^{xi)6{A-u) 



(23) 



Here the average is taken (i) over the positions of the maxima and of the envelopes ?/'^(x) and (2;), re- 
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spectively; (ii) over the statistical properties of the wave- 
function tails d/JxA[XAix)] and dfixBlXBix)] defined by 
Eq. ([7]), with the constraint Xajx) < [here we define 
Xaix) = ln|^Q(x)p, as in Eq. (|6j)]; (ni) over the energy 
difference e; and (iv) over the coefficient $ in Eq. (19 1. 

The following sections are devoted to extracting the 
three different regimes from the general formalism ( 22 ) 
and ( 23 ) : the behavior at the Mott scale and the leading 



and subleading corrections at sub-Mott lengths. 



V. BEHAVIOR AT THE MOTT SCALE 

As pointed out in the early work^^H^ ^jj^g hybridization 
of localized states introduces the logarithmically large 
"Mott scale" (|3| . The leading contribution to the behav- 
ior of R{uj,r) at r ~ Lm is obtained by picking out the 
following term from the general formula (l22|): 



iv+(xi)nv^_(x2)i^ 



iu+hV'A(xi)nv.B(x2)p. (24) 



Since the wave functions and ipB are localized at dis- 
tances of order one, and R{ijj,r) at r Lm varies at a 
logarithmically larger scale {5r ^ \/Lm, as shown be- 
low), the variables xa and xb are nearly pinned to the 
points xi and X2, respectively. Then integration over xa 
and Xb yields just the unit normalization of ipA and -0^, 
and we get 



= J dP{J) j d£5(A-cj)|u+|\ 



(25) 



where A and are functions of e and J defined by Eqs. 
(17) and (18). The measure of integration over J is 



dP{J) = f 



Po{xj,r)dxjdP{^), (26) 



where J is parameterized by Eq. (21), Po(xj,?') is the 
normal distribution (10 1, and /(xjAj is the cut-off func- 



tion [the same as in Eq. ([9|)]. The integral over e may be 
easily taken, which gives 



i?M(w,r) 



J 2 \uj e 



(27) 



where e = ^uP- — 4| J|2. 

Since the main contribution to the integral comes from 
logarithmically large intervals of xj, one can approxi- 
mate 

1 ~0(^_2|J|) (28) 



in Eq. (27) [the step function in the right-hand side takes 



care of the integration limits] and disregard the exact 
form of the distribution of <&. Then, within these ap- 
proximations, one gets 



i?A/(a;,r) 



2 In 



— Pm 
2^ 



1 + erf 



(29) 



i.e., the result reported in the last column of Tablejl] Note 
that the cut-off function fixj/f) does not play any role 
in this calculation, since /(—I) = 1 by the normalization 
of probability. 

We can further repeat the same procedure for the cor- 
relation function S'(a;,r) given by Eq. (23 1 by selecting 
the term 

'll}*^{xi)ll:+{x2)i^*^(x2)lp-[xi) 

-\u+\''\u^\^\4,a{x{)\^\^Pb{x2)? ■ (30) 

We then arrive at the formula similar to Eq. (25), but 
with replaced by — The formula (27) 

then gets replaced by 

SM{^,r)= j dP{J)\[^-'^) . (31) 

Now we cannot simply replace e by w, but need to expand 
to the next order. In fact, we can re-express 



1 / e uj 

2 Vcj ^ e 



d e 
dx.]^ 



-5[x.i-2\n 



-) 
2$/ 



which allows us to integrate over xj to obtain 



2^fwr 



exp 



4r 



(32) 



(33) 



again independently of the distribution of $ and therefore 
universally valid in SID and QID systems (including the 
numerical prefactoi'^. 



The result ( 29 ) has been prev iously rigorously der ived 
in SID and in QID - unitary case d^^ * '^^ , and the result ( 33 ) 
in the SID case^. Note that the location and width of 
the features in i?M('^j ^) and Snj{^, r) reflect directly the 
median and the width of the log-normal distribution for 
Xj in Eq. (21). In our ansatz (fTol), we take them related 



to each other, which cor responds to the single-parameter- 
scaling regim^^HmilllT] j^^ j^g£ ^ qualitative behav- 
ior of R{uj,r) at the Mott scale was also explained from 
the hybridization arguments, but the correct quantitative 
expression ( 29 ) could not be obtained without taking into 



account the log-normal distribution of the wave-function 
tails. 



VI. LEADING ORDER AT DISTANCES MUCH 
SHORTER THAN THE MOTT SCALE 

At distances 1 <C ?' <C Lm, the correlation functions 
R{uj,r) and S{uj,r) can be found, to the leading order, 
from the general expressions (22 1 and (23) if one retains 
only the contributions from ^a in both -0+ and 

|0+(a;i)|2|0„(a::2)|^ and V+(a;i)V'+(a;2)01(a;2)V'-(a^i) 
^ \u+\'\u.\'\i;A{xi)\'\i;A{x2)\' ■ (34) 



Then, using the relation (32 1, we can integrate over all 
the variables, except for xa and XA (in the order e, 
xb, $) and arrive at the result 

R{uj,r) w S{uj,r) 

2 / dxAdti.,AxA{x)]\i^Am^\'4'A{r)\^ ■ (35) 



Thus the short-distance behavior of these correlation 
functions is universal for SID and QID models and is 
only determined by the single-wave-function statistics. 
This result was rigorously derived for SID (as follows 
from Refs. [HI im and'^ and QlD-unitary casea^, and 
the exact form of this function is known. 



R{uj ^. 0, r) = Att' 



Qj.2 



tanhTrfc 

K CLrZ 7^ 6 

cosh nk 



(fe^ + l/4)r 



(36) 

Within our approximate method, we cannot derive this 
exact expression, but we can access its r 3> 1 limit. In 
this case, the main contribution comes from xa located 
between and r (see Fig. [3^) , with the two tails of the 
wave function ipA distributed independently, according 
to Eq. pll): 



r-O pO 

dxA / dxi I dx2 



X P{xi,XA)P{x2,r - XA)e^'^ 



-X2 



(37) 



If we use our ansatz ^ for P{x,r), then this integral 
formally diverges at xa ^ and xa — > r. This means 
that the main contribution to the correlation functions 
comes from configurations where the maximum of the 
wave functions coincides (within the localization length) 
with one of the two points. At such short distances, our 
ansatz (|9| is not applicable, but we can estimate the cor- 
relation function, up to a numerical prefactor, by cutting 
off the integral in xa within a localization length from 
and r [i.e., by integrating over xa in the limits {S, r — 5) 
with 6^1]. This immediately leads us to the asymptotic 
expression (at r 1) 



Riuj,r) « 5(w, r) cx r-3/2e-'-/^ 



(38) 



where the proportionality coefficient cannot be calculated 
within our approximation. 

The asymptotic ex pression (|38| ) is in agreement with 
the exact expression ( 36 IHHHI^ff Note that the form of 
the cut-off in the probability distribution ^ was impor- 
tant for calculating the correct power in the pre-exponent 
in Eq. ( 38 ). In fact, the correlation functions R{uj 0, r) 
0, r) are dominated not by "typical" local- 
but by the rare events, when the 



and S{lo 

ized wave functions, 
wave function ipA has two peaks at the positions and 
r of comparable height. This can also be seen from the 
exponential decay e~^^^, which does not describe the de- 
cay of a "typical" wave function [whose weight decays as 
e""", according to Eq. (12)], but is four times slower. 




FIG. 5: An illustration of notation for the calculation of the 
subleading order at 1 ^ r ^ Lm- The centers xa and xb 
of the pair of localized states are located outside the interval 
{xi, X2)- The variables ta, tb, and r used in Eq. (41 1 are the 
pairwise distances between the four points xa, x\, X2, and 
Xb- 



Somewhat similar rare events are important for the 
statistics of wave functions in the metallic limit^^, which 
also results in log-normal tails. However, the metallic 
regime is beyond the scope of the present paper. 



Note that the derivation of the relation ( 35 ) does not 
use the condition r ^ 1 (which is only needed for cal- 
culating its right-hand side), and is therefore valid for 
any r (in the limit w ^ 0). In terms of wave-function 



correlations, we may also rewrite Eq. (35) as 



R(uj 0, \xi - xal) = S{uj 0, \xi - X2\) 

= 2^1.-1 ^<5(S„ - E)\Mxi)\^\Mx2)\^ , (39) 

n 

(in this equation, we restore the physical units). This re- 
lation (without specifying the numerical prefactor) was 
already proposed in Ref. [33] based on similar hybridiza- 
tion arguments. However, the approach used in that work 
could not correctly reproduce the asymptotic behavior 
(38), since it did not include the log- normal distribution 
of the wave-function tails crucial for such a calculation. 



VII. SUBLEADING ORDER AT DISTANCES 
MUCH SHORTER THAN THE MOTT SCALE 

Remarkably, we can extend our method further to find- 
ing non-universal corrections SR{uj,r) and SS{Ltj,r) [de- 
fined in Eqs. (|4|] to the asymptotic behavior (38). Such 
corrections are given by the same "cross-terms" ( 24 ) and 
(30), as in the calculations of i?j\/(aj,r) and S'm(w7?') hi 
Section IVj One can check that, for this term, the main 



contribution comes from configurations with the points 
XA and Xb (the maxima of the wave functions ipA and 
^b) located outside the interval {xi,X2): see Fig. [5j 

For calculating 5R{u!,r), we start with the general ex- 
pression (22), where we only keep the term ( p4| ). Un- 
like in the calculation of Section |Vj here the tails of the 
wave functions contribute, and therefore we need to take 
into account their log-normal distributions. If we write 



= and \^Bix2)\' 



e^'^ , then the hy- 



bridization matrix element | J| may be expressed in the 
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form (21 1 with 



XJ = XA + XB+ X- 



(40) 



Here xAi a-nd xb are distributed with the probabiUty 
distribution ^ (with the cut-ofF) and the distribution of 
X is given by Eq. ( 10 ) (without the cut-off) [compare with 
an analogous form of Eq. ( 14 )] . In the calculations of this 



section, we further neglect the cutoffs, since all the three 
parameters — x, —XA, and ~xb are much larger than one 
and the integrals with respect to them can be done at the 
saddle-point level. The cut-ofF functions contribute only 
to the overall numerical coefficient, which is beyond the 
precision of our calculation. 

After integrating Eq. (22 1 over e, we arrive at 



SR{uj, r) (X drA / drs / dxA 

Jo Jo J -oo J-<x 

X PQixA,rA)Po{XB,rB)e^^+''" 
dxPo{x,r) I dP($)i(^ + 



dxi 



(41) 



Like in the calculation in Section [Vj we can use the ap- 
proximation ( 28 ). Furthermore, the integrals over ta and 



rg can be calculated at the saddle-point level (thereby 
fixing rA — —XA and — —Xb), and afterwards we 
integrate over XA and xb- The resulting expression is 



6R{uj,r) cx / dP($) 



dx 



exp 



4r 



X (z + l)e' 



2;=max(0,x-21n ^) 



(42) 



One can show that, at r < Lm/3, the main contribu- 
tion comes from the region x — 2 In ^ > 0. The inte- 
gral in X can be done in the saddle-point approximation, 
which sets x — ^3r. Finally, only the integral over $ 
remains [at our level of approximation, we also neglect 
-1-1 in the second line of Eq. (42 1]: 



2 2r 



6R{uj, r) cc (jj e 



dP($) 



In $ - In 



$2 



(43) 



where — exp [— (Lj\/ — 3r)/2]. Estimating the inte- 
gral (|43]) with the distributions dP{^) given by Eqs. ([20]) 
yields the results reported in the middle column of Ta- 
ble H 

An analogous calculation can also be performed for 



SS{iiJ,r) starting with Eq. (23 1. The calculation parallels 
the one above, with the onl; 



y difference that {l/2){e/u}- 



uj/e) in Eq. (41) must be replaced by 

1 / e uj 
2\uj e 
This results in 



2 In 



2$ 



X-Xa-Xb 



(44) 



5S{uj,r) cx - / dP($) 



dx 



exp 



ix + ry 

4r 



X ze 



(45) 



To the precision of our approximation, this expression 
is opposite in sign to Eq. (42). Therefore, we conclude 
that, within our approximation, 5S{ui^r) « —SR{uj,r). 
The corresponding formulas are reported in the middle 
column of Table [llj Note that our method does not give 
the numerical coefficients in dR{Lu,r) and SS{uj,r), but 
predicts that they have the same absolute value and are 
opposite in sign [positive for dR{Lu,r) and negative for 
SSiio,r)]. 

One can compare our results of this section with the 
exact calculations. The only case, where a direct compar- 
ison is possible is the QlD-unitary case, where SR{uj,r) 
was computed in Ref. iSOj and, to the leading order, co- 
incides (up to a numerical prefactor) with our present 
result. Note that our results for 5R{uj,r) and SS{uj,r) 
in the SID and QlD-uintary cases are only applicable at 
r < Lj\//3. The new length Lm/3 appeared in Ref. [30] 
as the distance at which the correction 6R(uj,r) in the 
QlD-uintary case changes its asymptotic form (in tech- 
nical terms, there was a switching of the pole and the 
saddle in the integral determining the leading form of 
the correction). 

Another situation where an indirect comparison can 
be made is the SID case. There, using the formalism of 
Ref. [27l one can sho'wi^that 5R{uj,r), at small r starts 
with the order Inw, i.e. consistent with our result for 
the SID case (strictly speaking, the comparison is not 
accurate, since our result is only applicable at r ^ 1 
while the expansion of the result from Ref. [27| is done 
at r <C 1 but we expect that the leading w dependence 
is the same in both regimes). In the other cases, there 
are no exact calculations to which our results reported in 
Tables |l] and |ll] could be compared, thus they should be 
considered as conjectures. 



VIII. SUMMARY, DISCUSSION, AND 
OUTLOOK 

To summarize, in this paper we propose a simple tech- 
nique of treating correlation of wave functions in An- 
derson localization in terms of hybridization of local- 
ized states. While the essence of our method repeats 
the well-known Mott argument, supplementing it with a 
log-normal probability distribution for the tails of local- 
ized wave functions (and, consequently, for the hybridiza- 
tion matrix elements) gives the method a quantitative 
power. We have checked that the results produced with 
our simplified method reproduce quantitatively the main 
features of the available exact results (obtained by more 
sophisticated techniques). 

Nevertheless, we should emphasize that the presented 
method remains a phenomenological recipe, and its justi- 
fication still needs to be completed. The method depends 
on several assumptions of various level of rigor. For the 
benefit of the reader, we list them below: 

1. A possibility to define localized wave functions that 



are further hybridized into the eigenstates ( 16 ) 
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close in energy. These wave functions tpA and ■0^ 
are not rigorously defined in our argument, and the 
formalization of this step would be helpful for a rig- 
orous justification of the method. 

2. The log-normal distribution of the wave-function 
tails (10), supplemented by a suitable cutoff (|9|. 
While we present some arguments in favor of these 
formulas, they are not formally derived. We hope 
that a rigorous derivation of this step may be pos- 
sible with the methods of Ref. 9, 



3. The hybridization matrix elements (21) are as- 



sumed to be proportional to i})a{x)'iI)b{x)^ the prod- 
uct of the two hybridizing tails, and therefore to 
obey the same log-normal distribution. Since nei- 
ther J nor ^ja.b are formally defined, this assump- 
tion also remains a phenomenological construction. 



4. The factor $ in (21) reflecting the interference of 



channels. Its probability distributions (20 1 are in- 
troduced phcnomenologically. 

Note that, for different calculations, different assump- 
tions play a role. At the Mott length scale (Section [V]), 
we only used the assumptions 1 and 3, with the assump- 
tions 2 and 4 being irrelevant. For the leading behavior 
at \ <^ r <^ Lm (Secti on |VI[ ), we also used the assump- 
tion 2, while in Section |VII| we additionally need the as- 
sumption 4 to calculate the subleading terms SR{uj,r) 
and SS{u!, r). 

We wish to remark here that our method was devel- 
oped under the specific assumption of a weak disorder 
(specifically, it is applicable in the model of the Gaussian 
white noise), which implies the single-parameter-scaling 
relation (12 1 between the variance and the average of 
the wave-function tail. However, the method can be ex- 
tended to other interesting models of localization by re- 
laxing this assumption. This would imply a modification 
of the log-normal probability distributions ^ and ( 10 ) 
in the assumptions 2 and 3. One example where such a 
modification may be applicable is localization far below 
the mobility edge (see, e.g., Refs. l40H42j) . Another inter- 
esting example is the exactly solvable localization prob- 
lem with the Cauchy-distributed disorder considered in 
Ref. 26, In that model, the single-parameter-scaling re- 
lation (12) does apply, but with a different coefhcient. 



Finally, it would be interesting to extend our approach 
to higher dimensions. The key issue for such an exten- 
sion is the log-normal distribution of the wave-function 
tails (our assumption 2 above). While we are not aware 
of such results for wave functions in higher-dimensions, 
a similar claim was made for the probability distribu- 
tion of the conductance^. Namely, it was shown that, in 
the weak-scattering case and in the insulating phase, the 
conductance distribution in the large-system-size limit 
is log-normal with the single-parameter-scaling relation 
between the average and the variance of the form (12 1. 



We believe that such a system can also be treated by 
our method with a suitably modified probability distri- 
butions ^ and (flO|. 



Therefore one may assume that a similar universal dis- 
tribution is also valid for the wave-function tails. If it is 
indeed the case, then our calculations in Sections |V] and 
|VI| can be straightforwardly extended to higher dimen- 
sions. In particular, the counterpart of Eq. (39) in any 
dimension would read 

R{uj 0, \xi - X2I) = S{uj 0, l^i - X2\) 
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(46) 

where Sd-i is the area of the [d— l)-dimensional sphere 
(e.g., in our one-dimensional case, = 2), and the 
definition ([3]) of Lm remains the same in any dimen- 
sion. Also, under the same conditions, the behavior of 
_Ra/ (i^, 1") and Sm{i^, r) (the right column of Tables|l]and 
|ll] calculated in Section |v]) would be universal in any di- 
mension. This implies, in particular, the validity of the 
Mott formula for the frequency-dependent conductivity 
(j{uj) oc a;^(lnti;)''+^ in any dimension (which can be de- 
duced from S'M(w,f), see, e.g., Ref. f28|). in agreement 
with the original argumen'p2l. 
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